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Abstract 

It is shown that existing data on the mixing between up and down 
fermion states and on the hierarchical mass ratios between fermion 
generations, as far as can be so analysed at present, are all consistent 
with the two phenomena being both consequences of a mass matrix 
rotating in generation space with changing energy scale. As a result, 
the rotating mass matrix can be traced over some 14 orders of magni- 
tude in energy from the mass scale of the i-quark at 175 GeV to below 
that of the atmospheric neutrino at 0.05 eV. 



1 Introduction 



Along with the mystery of why there should be in nature 3, and apparently 
only 3, generations of fermions, the fact that their masses should be hierar- 
chical and that they should mix, as embodied for quarks in the CKM matrix 
IJ and exhibited for leptons in neutrino oscillations, has remained one of 
the great puzzles of particle physics. In the context of what we call the 
Dualized Standard Model (DSM) which is an explicit attempt to solve 
this generation puzzle, we have suggested that the mass hierarchy and the 
mixing phenomenon can both result from a mass matrix which changes its 
orientation in generation space (rotates) with changing energy scale and have 
obtained rather good agreement with experiment based on this hypothesis. 
This previous treatment, however, does not fully cover the following two as- 
pects. Firstly it depends on the details of the DSM mechanism driving the 
rotation, which details may not be strictly necessary for deriving the said 
result, and secondly, it does not clearly reveal the degree of significance of 
the claimed agreement with experiment nor the amount of direct empirical 
support, if any, for mass matrix rotation. For this reason, our purpose in 
this paper is to leave aside for the moment the driving mechanism (whether 
DSM or otherwise) but put instead the rotating mass matrix directly to the 
test by going straight to the experimental data and seek evidence there for 
the rotation hypothesis. 

We shall show that by inputting all the available mass and mixing data 
on both quarks and leptons, assuming only that these all arise from mass 
matrix rotation in the manner to be explained below, one can trace the 
implied rotation over a scale range of some 14 orders of magnitude. The 
result is seen to be fully consistent with the fermion states all lying on a 
single smooth rotation curve linking the t quark at 175 GeV through all the 
intermediate fermion states down to the second heaviest neutrino vi at less 
than 10~ 2 eV. This constitutes positive evidence for mass matrix rotation, 
although as yet only circumstantial, which is independent of the theoretical 
mechanism responsible for driving it, whether DSM or otherwise. It will be 
seen, nevertheless, that the shape of the empirical rotation curve so traced is 
indicative of having the two rotational fixed points at infinite and zero energy 
scales predicted by and, as far as we know, special to the DSM scheme. 

We begin by re-examining the reasoning behind the ansatz that both mass 
hierarchy and mixing can arise from a rotating mass matrix, for although it 
has been stated here and there and used already in several earlier papers, we 
cannot assume that it is of general knowledge. The first point to note is that 
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once the mass matrix rotates with changing scale, then even such familiar 
concepts as particle masses, state vectors and mixing of flavour states, will 
have to be revised or refined, as we shall now explain. 

Let us start with a fermion mass matrix traditionally defined by a term 
in the action of the form: 

M + t-C, (1) 

where ip L and ip° R represent respectively the left- and right-handed fermion 
field, each being a vector in 3-dimensional flavour space, here given in the 
weak gauge basis, and m is a 3 x 3 (complex) matrix. The matrix m can 
always be diagonalized thus: 

U\mU R = diag{mi, m 2 , m 3 } (2) 

with Ul, Ur unitary and m, : real. Thus in terms of the fields: 

$L = uhl>k ^ R = uy R , (3) 

the term (H) in the action takes on the diagonal form: 

^ L diag{mi, m 2 ,m 3 }i/; R . (4) 

When the mass matrix m is constant in orientation with respect to scale 
change, i.e. in the language of this paper, when the mass matrix does not 
rotate, which is the simple case usually considered, then the particle masses 
of the 3 flavour states are just given by the diagonal values m^. The above 
apply to both up and down quarks in the case of quarks, and to both charged 
leptons and neutrinos in the case of leptons. Hence, from the mass matrix, 
one obtains for the up and down states each a diagonalizing matrix Ul which 
we can denote respectively as Ul and U' L . Again, in the simple case when the 
mass matrices do not rotate, then the mixing matrix between up and down 
states (i.e CKM jl] for quarks and MNS H for leptons) is given just as ||]: 

V = U L U'l (5) 

For the discussion in this paper, it is more convenient to work with an 
equivalent form of the mass matrix adopted by Weinberg in ||. Since the 
right-handed fermion fields are flavour singlets, they can be arbitrarily rela- 
belled without changing any of the physics. (Witness the fact that the mixing 
matrices between up and down states depend only oyiUl not on Ur.) Hence, 
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by an appropriate relabelling of right-handed fields, explicitly by defining 
new right-handed fields: 

■K = u L uy% (6) 

one obtains ([!]) in a form in which the the mass matrix becomes Hermitian: 
$m w ]-{l + 7 5 )V> + $m w ~(l - 75)^ = 4>m w ij, (7) 

with 

mw = mUjiU^. (8) 

This is convenient because in the simple case when the mass matrix does not 
rotate, the particle masses are now just the real eigenvalues of the Hermitian 
matrix mw, as can readily be checked with (||). Furthermore, the mixing 
matrix between up and down states become just 

V ij = {v i \y f j ), (9) 

with |vj) being the eigenvector of mw corresponding to the eigenvalue mi, 
and the prime denoting down-type fermions. The scalar product (Vj|v'j) is 
of course an invariant independent of the frame in which these vectors are 
expressed. Thus, in terms of this Weinberg form of the mass matrix mw 
(which in this paper will be used exclusively and from which the subscript 
W for convenience will henceforth be omitted), the mass values and state 
vectors of flavour states, as well as the mixing matrix between up and down 
fermions, are all easily defined in the case of no rotation in the mass matrix. 

Consider now what happens in the case when the mass matrix rotates 
with changing scale as examined in this paper. Both its eigenvalues and their 
corresponding eigenvectors now change with the scale so that the previous 
definition of these as respectively the masses and state vectors of flavour 
states is no longer sufficiently precise, for it will have to be specified at which 
scale(s) the eigenvalues and eigenvectors are to be evaluated. 

In the simple case of a single generation, i.e. when the the mass matrix is 
just a number, one is used to defining the particle mass as the running mass 
taken at the scale equal to the mass value itself, i.e. at that /x when \i = m(fi). 
One might be tempted therefore to suggest for the multi-generation case that 
one defines the mass m; and the state vector of the state i, as respectively 
just the ith eigenvalue and eigenvector of the matrix m taken at the scale 
Hi = mj(/ij), with mj(/i) being the scale-dependent ith eigenvalue of the 
matrix m. However, such a definition will not do, because it would mean that 
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the state vectors for the different generations i will be defined as eigenvectors 
of the matrix m at different scales. Although the eigenvectors % for different 
eigenvalues % are orthogonal, m being Hermition, when taken all at the same 
scale, they need not be mutually orthogonal when taken each at a different 
scale. But the state vectors for different flavour states ought to be orthogonal 
to one another if they are to be independent quantum states. Otherwise, it 
would mean physically that the flavour states (such as the charged leptons 
e, n and r, for example) would have nonzero components in each other and be 
thus freely convertible into one another, or that the mixing matrices (CKM 
for quarks and MNS for leptons) would no longer be unitary, which would of 
course be unphysical. 

How then should the mass values and state vectors of flavour states be 
defined in the scenario when the mass matrix rotates? To see how this 
question may be resolved, let us examine it anew with first the [/-type quarks 
as example. The 3x3 mass matrix m has 3 eigenvalues with the highest 
value rrt\ corresponding to the eigenvector Vi, both depending on scale 11. 
Starting from a high scale and proceeding lower, one reaches at some stage 
\i\ = mi(//i), i.e. when the scale equals the highest eigenvalue m\. One 
can then naturally define this value m 1 ( / u 1 ) as the t quark mass m t and the 
corresponding eigenvector Vi(fii) as the t state vector v t . Next, how should 
one define the mass m c and the state vector v c ? We have already seen above 
that they cannot be defined as respectively the second highest eigenvalue rri2 
of the 3x3 mass matrix m and its corresponding eigenvector at the scale 
/i2 = ^2(^2), because this vector is in general not orthogonal to the state 
vector v t which the state vector v c ought to be. It is not difficult, however, to 
see what is amiss. At scales below the t mass, i.e. when \i < m t , t would no 
longer exist as a physical state, so that what functions there as the fermion 
mass matrix is not the 3x3 matrix m but only the 2x2 submatrix, say rh, 
of m in the subspace orthogonal to v t . Hence, for consistency, one should 
define m c as the highest eigenvalue rfi2 of the matrix rh and the state vector 
v c as the corresponding eigenvector, both at the scale £12 = ^(p^)- The 
state vector of c so obtained is of course automatically orthogonal to v t as 
it should be. Repeating the argument, one defines further the mass m u and 
state vector v u respectively as the "eigenvalue" and "eigenvector" of rh at 
the scale fi 3 = rh 3 (fi 3 ), with rh being the 1 x I submatrix of m in the subspace 
orthogonal to both v t and v c . Proceeding in this way, all masses and state 
vectors are defined at their own proper mass scale and the state vectors are 
mutually orthogonal as they should be. Besides, though stated above only 
for 3, the definition can be extended to any number of fermion generations, 
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should there be physical incentive for doing so. 

We note that a refined definition of masses and state vectors for flavour 
states is a question which has to be adressed in principle when the mass ma- 
trix rotates with changing scale, whatever the speed of the rotation. Hence, 
it ought to figure even in the conventional formulation of the Standard Model 
where the mass matrix is bound to rotate when there is nontrivial mixing 
between up and down states || fjj , although the rotation there is rather slow 
and its effects are for most practical applications negligible. The above solu- 
tion to the problem, though first made in the context of the DSM scheme f7|, 
is seen actually to apply to any rotating mass matrix, and it is, as far as we 
know, the only solution to the question yet given in the literature. So this is 
the definition we shall employ in what follows to analyse existing mass and 
mixing data in terms of a rotating mass matrix. The success or otherwise of 
the analysis could thus itself be regarded as a check on the validity of the 
above solution. 

Having now made clear our general procedure for defining masses and 
state vectors, let us return to the problem at hand, namely that of testing 
empirically the hypothesis that both mass hierarchy and mixing arise as 
consequences of the rotation of the mass matrix. First, by 'mixing arising 
from rotation', we mean that the mass matrices of up and down fermions 
taken at the same scale are aligned in orientation at all scales, only differing 
by a normalization factor, and it is the rotation alone which is giving rise to 
the mixing. Secondly, by 'mass mass hierarchy arising from rotation' as well, 
we mean that the mass matrix has at all scales only one massive eigenstate, 
with the masses of the lower generations appearing only by virtue of the 
rotation via what we called the "leakage" mechanism to be explained below. 
We have thus for both up and down fermions a mass matrix of the form: 



where my is the single nonvanishing eigenvalue of m and |r) its corresponding 
(normalized) eigenvector. Or explicitly for 3 generations, we have: 



where only m,T depends on the fermion species. Thus the whole content 
of the rotating mass matrix is now encapsulated in r = (£, 77, C)> a vector 
rotating in 3-dimensional generation space with changing scale \x. 



m = mr|r) (r 



(10) 




(11) 
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That such a simple form of the mass matrix can nevertheless give rise 
to nontrivial mixing and mass hierarchy is most easily seen in a simplified 
situation when there are only 2 generations instead of 3. When applied to 
this simple case, the procedure detailed above for defining masses and state 
vectors of flavour states gives the state vector of t as the single massive 
eigenstate r of the [/-quark mass matrix at the scale \x = m t , and the vector 
v c as a vector orthogonal to v t , as depicted in Figure [l]. Thus v c has by fllOD 
a zero eigenvalue for m at the scale \i = m t . But according to our previous 
conclusion, this should not be interpreted to mean that c has a zero mass. 
The c mass m c is given instead as the eigenvalue of m at scale /i = m c , 
which in this simplified 2 generation case is just the expectation value of m 
in the state v c . At the scale fi = m c , however, the vector r is already rotated 
to a direction different from that of v t , as illustrated in Figure [l], and has 
acquired a component equal to sin^c in the direction of v c , with 9 tc being 
the rotation angle between the scales [i = m t and = m c . Hence, according 
to the above definition, one obtains "by leakage" a nonzero mass for c given 



by the expectation value of ( JT0| ) with respect to v c , namely: 

m c = m t sin 2 9 tc . (12) 

Similarly, although the mass matrices of the U and D quarks according to 
( |T0| ) are aligned in orientation at all scales, one sees from Figure [| that by 
virtue of the rotation of the vector r from the scale \x = m t to the scale 
fj, = m b where the state vectors v t and v b are respectively defined, the two 
state vectors will not point in the same direction. One has thus from (^) 
simply by virtue of the rotation a nonzero mixing between the t and b states 
with the CKM matrix element: 



V t b = v t .v 6 = cosOtb 7^ 1, (13) 

where 9 tb is the rotation angle between the two scales. 

One sees already from these examples in the simplified scenario of only 

2 generations that both lower generation masses and nontrivial mixing will 
automatically be obtained from a rotating mass matrix even if one starts 
with neither. One sees also that the same conclusion can be drawn in the 

3 generation case so long as the masses and state vectors of flavour states 
are defined for a rotating mass matrix as detailed above. The remaining 
question to ask is then whether the actual masses and mixing of fermion 
states as observed in experiment can actually be explained in this simple 
way from a rotating mass matrix, and this is the question we set out to 
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Figure 1: Masses for lower generation fermions from a rotating mass matrix 
via the "leakage" mechanism. 

answer in this paper. This question can be adressed empirically since, as 
we shall show, many of the relevant quantities have already been measured 
and need only to be arranged and interpreted appropriately for the present 
purpose according to the definitions detailed above. 

2 Analysis in the planar approximation 

Let us first perform the analysis in the simplified situation with only the 
2 heaviest generations in each fermion species, which scenario makes the 
analysis much more transparent since the problem then becomes planar and 
there is only one real rotation angle to consider |4[]. This simplification will 
be shown later to approximate already rather well the actual 3-generation 
situation in the scale region above roughly the [i lepton mass. We have then 
the pictures shown in Figures [I] and ^] for obtaining respectively the lower 
generation masses and mixing matrix elements. 

Consider first mixing matrix elements. Suppose from the scale of the t 
mass to that of the b mass, the mass matrix has rotated by an angle 6tb 
(Figure [|). As explained above, one easily obtains then the CKM elements 
as: V t b = cos 9 t b and \V ts \ = \V c b\ = sin 6^. From the measured values of these 
elements given in the latest databook ||, namely: 

| Kb I = 0.9990 - 0.9993, \V ts \ = 0.035 - 0.043, \V cb \ = 0.037 - 0.043, (14) 
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Figure 2: Mixing between up and down fermions from a rotating mass matrix. 

one gets thus from each an estimate of the rotation angle, respectively: 

6a = 0.0374 - 0.0447, 0.0350 - 0.0430, 0.0370 - 0.0430, (15) 

the values obtained being fully consistent with one another. (One notes that 
from the same Figure 0, one could deduce in principle also V cs = cos 9 t b, 
but this will be seen, in contrast to the 3 other mixing elements already 
considered, to be a poor approximation receiving large nonplanar corrections 
when all 3 generations are taken into account.) 

Consider next the second generation masses obtained by the leakage 
mechanism. Suppose from the scale of the t mass to that of the c mass, 
the mass matrix has rotated by an angle 8 tc - Using ( |T2"D and the measured 
values of m t and m c given in ||, namely: 

m t = 174.3 ±5.1 GeV, m c = 1.15 - 1.35 GeV, (16) 

one obtains the estimate: 

tc = 0.0801 - 0.0894. (17) 
Similarly, from the measured values from ||: 

m b = 4.0 - 4.4 GeV, m s = 75 - 170 MeV, (18) 
one obtains the estimate: 

6s = 0.1309 - 0.2076, (19) 
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the error being so large because of the intrinsic uncertainty in defining the s 
quark mass, while from the measured values from ||: 

m T = 1.777 GeV, m M = 105.66 MeV, (20) 

one obtains the estimate: 

d Til = 0.2463. (21) 

Assume now that the mass matrices of the U and D quarks as well as 
the charged leptons are all aligned at the same scale as proposed above, 
and plot the values of the rotation angles obtained before, all starting from 
the direction of the t quark state. One obtains then the Figure |3] where, 
rotation angles being simply additive in the planar approximation, we have 
taken 6 ts = 9 t b + 9b s , and 9 tfJL = 9 tr + 9 Tfl , with 9 t b taken from flip] ) and 
9 tT (indicated by a cross in Figure |3|) estimated by interpolation between the 
values of 6 tb and 9 tc given above. One sees first that the information gathered 
so far indeed appear consistent with the data points lying on a single smooth 
rotation curve as can be explained by a mass matrix rotating smoothly as 
the scale changes, which we regard as already a nontrivial support for the 
rotation hypothesis independent of any other theoretical consideration. In 
addition, one notes that the data by themselves are already indicative of the 
rotation angle approaching an asymptotic value thus suggesting a rotational 
fixed point at infinite scale as predicted by DSM. Indeed the data points all 
sit surprisingly well on the rotation curve obtained a few years earlier |J in 
our DSM calculation. Making a best fit to the data by MINUIT produces 
an excellent fit: 

9 = exp(-2.267 - 0.509 ln/i) - 0.0075, (22) 

for fi in GeV, with a \ 2 °f 0.21 per degree of freedom, which is hardly 
distinguishable from the DSM curve, as can be seen in the same figure. 

The above analysis was performed under the simplifying assumption of 
there being only 2 generations of fermion states but one can show that for the 
quantities so far considered it is already a good approximation to the actual 
3-generation situation. When all 3 generations are taken into account, the 
mixing matrix can be parametrized as: 

' V tb V ts V td \ ( ci -sic 3 -S1S3 \ 

V ch V cs V cd = sic 2 C1C2C3 - s 2 s 3 e lS C1C2S3 + s 2 c 3 e lS , (23) 
V V uh V us V ud ) \ S1S2 C!S 2 c 3 + c 2 s 3 e lS C1S2S3 - c 2 c 3 e lS ) 
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Figure 3: The rotation angle changing with scale as extracted from data on 
mass ratios and mixing angles and compared with the best fit to the data 
(dashed curve) and the earlier calculation by DSM (full curve) 0, in the 
planar approximation. 
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which is the original Kobayashi-Maskawa parametrization ]IJ], only with the 
ordering of fermion states reversed. One sees then that if we continue to 
denote as before V t b as cos 6^, the elements V ts and V c b are no longer just 
given by sm9 t b but by respectively sin Q t b cos #3 and sin 6 t b cos 62- The angles 
62 and #3, however, are small as can be estimated from the empirical values 
given in || for the corner elements of the CKM matrix in comparison to the 
values of V ts and V c b quoted above, giving: 

= 0.004- 0.014 — > I tan 9 3 \ = 0.093 -0.400, 
I Kb I = 0.002- 0.005 — ► |taa0 2 | = 0.047-0.135, (24) 

from which one gets: 

cos 9 2 = 0.999 - 0.991; cos 6 3 = 0.996 - 0.928. (25) 

Hence, one concludes that in the 2 generation planar approximation of Fig- 
ures and ^ where one puts V ts = V c b = sxn9 tb , one has made an error of 
at most a few percent which is seen hardly to affect the plot or any of the 
previous remarks we made. A similar error has been made in Figure |3] as 
regards the lower generation masses obtained from the "leakage" mechanism. 
The estimate ( |i~9"D is for the angle rotated between the scales of mj and m s 
but in Figure |] we have added this angle to the rotation angle from scale 
m t to scale m& to get the angle from scale m t to scale m s . Such an addition 
is valid in the 2 generation approximation but has nonplanar corrections in 
the actual 3 generation situation. The error so incurred can be estimated as 
follows. The angle between the plane defined by the t and c vectors and the 
plane defined by b and s is given by the angle between their normals, namely 
the vectors for u and d respectively, which according to || takes the value: 

\V ud \ = cos6 ud = 0.9742 - 0.9757 (26) 

giving 

ud = 0.2209 - 0.2276. (27) 

The nonplanar error incurred in the angle at scale m s plotted in Figure |3| is of 
order 1 — cos 9 u d and is thus of order again a few percent, which is negligible 
given the large error already inherent in the definition of the s quark mass. 
A similar error is presumably present in the angle plotted in Figure ^| at scale 
m /t . More details of this analysis can be found in a preliminary report ||10| . 

However, though good for illustration purposes because of its simplicity, 
the above analysis in the planar approximation is incomplete in that it is 
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restricted only to scales above the /x mass where nonplanar effects are neg- 
ligible, and hence cannot account for all the available empirical information 
which includes the data from neutrino oscillations at very low scales where 
nonplanar effects can no longer be neglected. Furthermore, to test the rota- 
tion hypothesis exhaustively one has also to ensure that no hidden violation 
of the hypothesis exists in the off-planar direction. In any case, the planar 
approximation, though good for illustrative purposes, actually need not be 
made since, as detailed in ([H]) above, the whole content of the rotating mass 
matrix even when all 3 generations are taken into account is encapsulated 
just in the single 3-dimensional vector r(/x) depending on scale \x. The only 
technical problem posed by a full analysis is thus how to extract this vector 
at various scales from the existing data on fermion mass ratios and mixing 
parameters. Once so extracted, the vector can be confronted with the rota- 
tion hypothesis and should be consistent with tracing out a continuous curve 
in 3-space if the hypothesis is correct. For the extraction of this vector, we 
now proceed as follows. 



3 Extracting r(/i) from quark data 

The ?7-quarks t, c, u are independent quantum states so that their state vec- 
tors should form an orthonormal triad in generation space, which we can 
choose without loss of generality as: 

v t = (1,0,0); v c = (0,1,0); v„ = (0,0,1). (28) 

The .D-quark state vectors also form an orthonormal triad the orientation of 
which relative to the [/-triad is given by the CKM matrix elements: 

v 6 = (V tb , V cb) V ub ); v s = (V ts , V CS1 V us ); v d = (V td , V cd) V ud ). (29) 

Hence, if the complex elements of the CKM matrix are accurately known, the 
.D-triad would also be determined. At present, however, only the absolute 
values of the CKM matrix elements are experimentally known to reasonable 
accuracy, leading thus to some ambiguities in the determination of the D- 
triad. In particular, one is forced to ignore for the moment in the CKM ma- 
trix the CP-violating phase which is experimentally very poorly determined 
and treat the D-triad also as real vectors. Inserting then the experimental 
limits on the CKM matrix elements as read from || gives rather tight con- 
straints on the directions of the D-triad with errors so small as to be mostly 
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negligible for our analysis: 



v fe = ((0.9990 - 0.9993), (0.037 - 0.043), -(0.002 - 0.005)); 
v s = (-(0.035 - 0.043), (0.9734 - 0.9749), (0.219 - 0.226)); 
v d = ((0.004 - 0.014), -(0.219 -0.225), (0.9742 - 0.9757)). (30) 

The signs of the 3 components for b can be chosen arbitrarily by choosing the 
physically irrelevant phases of the various quark state vectors and a particular 
choice has been made in (|3"0| ) for convenience. The signs for the other 2 states 
are then determined by orthogonality. Actually, given the present errors on 
the CKM matrix elements, there is an alternative solution to that shown for 
the state vectors of s and d, which however we can ignore for reasons to be 
explained below. 

We notice that by the considerations in the preceding section, the state 
vector of b, this being the heaviest state in the D sector, is just the rotating 
vector r(/i) taken at the scale fM — rrib, thus: 

r(m&) = v 6 . (31) 

Together with: 

r(m t ) =v t = (1,0,0), (32) 

we have then two points on the trajectory for r(/x) we wish to trace. One 
convenient way to present this, we find, is to write this rotating vector as 
r 0") = (£(/■*), ^(/^CCaO) with £(/^) 2 = 1 - "nW - C(a0 2 and plot the con- 



straints on t}(ji) and on the r/^-plane. The result from ( p2|) and (|3lD 
are then entered as the first 2 points from the left in the 3-D plot of Figure 
U This plot, which shows rj(fi) and as functions of the energy scale \x 
can in principle incorporate all the information that we shall extract from 
data. However, it being often hard to read the information it contains, we 
shall supplement it by its 3 projections onto the 3 co-ordinate planes, namely 
onto the ?y£-plane in Figure H the /x?7-plane in Figure ||, and the /z£-plane 
in Figure [7|, which projections, as we shall see, will be useful later also for 
interpolation and extrapolation purposes. 

Next, we turn to consider r(/x) at the scale \i = m c . Since the mass 
of m c , by our original hypothesis, comes about only through the "leakage" 
from m t as detailed above, it follows that the vector r(m c ) will have to be 
a linear combination: cos 9tcVt + sin^ c v c , of the state vectors v t and v c of 
respectively t and c with sin# te given by (]I2]), or explicitly: 



'(m c ) = \Jl- m c /m t v t + \J m c j m t v c . (33) 
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Figure 4: A plot of the rotating vector r(/x) as extracted from existing data 
on fermion mass ratios and mixing parameters, where its second and third 
components, i.e. rj(fi) and C^); are plotted as functions of ln/i, /i being 
the energy scale. The experimentally allowed values at any one scale are 
represented as an allowed region on a plaquette, with the scale corresponding 
to a plaquette being given by the intersection, denoted by a small circle, of its 
left-most boundary with the /i-axis. For example, the first small plaquette on 
the left of the figure corresponds to the scale /i = m b , on which plaquette the 
allowed region for r(/i) = v?, is very small because of the small experimental 
error on the CKM matrix elements Vtb, V c b and V u b. The last plaquette on 
the right, on the other hand, corresponds to the scale fi = m Ua , on which 
plaquette the allowed region for r(/i) is a rough rectangular area bounded by 
the data on v oscillations from atmospheric neutrinos and from the Chooz 
experiment. The curve represents the result of a DSM one- loop calculation 
from an earlier paper |J which is seen to pass through the allowed region 
on every plaquette except that for the electron e. For further explanation of 
details, please see text. 
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Figure 5: Projection of Figure f| onto the ^-plane. The full curve repre- 
sents the DSM one-loop calculation of |J and the dashed curve its suggested 
deformation at low scales to fit the data on m e and U e 2- 
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Figure 6: Projection of Figure |] onto the /x^-plane. The full curve repre- 
sents the DSM one-loop calculation of |J and the dashed curve its suggested 
deformation at low scales to fit the data on m e and U e 2- 



17 




Figure 7: Projection of Figure |] onto the /i£-plane. The full curve repre- 
sents the DSM one-loop calculation of |J and the dashed curve its suggested 
deformation at low scales to fit the data on m e and U e 2- 
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Inputting the quoted experimental limits on m t and m c gives us the allowed 
region for the third data point on the trajectory of r(/i), as is shown on 
Figure There is in fact of course another solution for r(m c ) in ( |3"3"D with 

sm9 tc = —\Jm c /m tl meaning that there is an additional disjoint branch of 
the allowed region in another quadrant of the 7^-plane to that displayed 
in the figure. The existence of this additional branch, however, does not 
affect the question of interest to us here, namely, whether the allowed region 
is consistent with data lying on a smooth rotation curve, so long as the 
first branch already does, and can therefore be ignored. Such additional 
branches of the allowed regions will in fact always occur for all other pieces 
of information on r(/x) that we shall extract from data, since experiment so 
far gives only the absolute values of the relevant quantities but, for the same 
reason as above, these additional disjoint branches for our present purpose 
can almost always be ignored. 

The next point in line is r(/i) at /i = m s , which will of necessity be poorly 
determined because the s mass is very poorly known. Nevertheless, whatever 
is taken for the mass of s so long as it is given by the "leakage mechanism" 
from b, r(m s ) will have to be a linear combination of the state vector v b of b 
and the state vector v s of s, so that by reasoning exactly as above for c, we 
have: 

r(m s ) = \Jl- m s /m b v b + \Jm s /m b v s . (34) 

The range of values for m s is given in || as 75 — 170 MeV, and in |ll| as 
100 — 300 MeV, depending on the scale at which the limits were determined 
respectively. The allowed region shown for r(m s ) in Figure |] corresponds 
to the union of the above two ranges for m s values, which allowed region, 
one will notice, does not lie entirely on the plaquette corresponding to an m s 
value at 176 MeV but protrudes to either side of it. This gives then a 4th 
point on the trajectory for r(/x). As noted before, given the present errors on 
the CKM mattrix elements, there is actually an alternative solution for v s 
with a different sign for the third component to that given in (|3T!D, meaning 
an additional branch to the allowed region for r(m s ), but this can be ignored 
for the same reason as that given for r(m c ) above. 

The above represents more or less the most that can be extracted from the 
data on quark masses and mixing (barring CP violation) about the rotating 
vector r(n) apart from some as yet rather uncertain information from the 
masses of the light quarks u and d. The method we used for extracting the 
fermion masses from the "leakage mechanism" as contained in e.g. (|12|) was 
meant only for freely propagating particles and should not in principle be 
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applied to quarks which are confined, except approximately to the heavier 
quarks which are generally regarded as quasi- free. For the light quarks u and 
d which are tightly confined, it is clearly not applicable. The question then 
arises in what way these light quark masses are to be defined. Experimentally, 
these masses are determined at some convenient but somewhat arbitrary scale 
such as 1 or 2 GeV, and it is not clear what these values should correspond to 
in the "leakage mechanism" . One possibility is to consider these mass values 
as the "leakage" from the rotating vector r(/x) taken at the chosen scale 1 or 
2 GeV into the vectors v u and v rf thus: 



In that case one obtains values for both m u and of the right order of 
magnitude in the MeV region, but it is not certain whether this is of much 
significance. As matters stand, therefore, we can only leave open the question 
of the light quark masses. 

4 Extracting r(/i) from lepton data 

To proceed further, one turns now to the leptons for which the preceding 
analysis for quarks can in principle be independently repeated, for as far as 
the premises of the rotating mass matrix as set up at the beginning of this 
paper is concerned, there is strictly nothing which needs connect the mass 
matrices of quarks and leptons. However, the DSM scheme suggests that the 
rotating matrices for quarks and leptons both lie on the same trajectory, the 
last being specified just by the vev's of the (dual colour) Higgs fields which 
are independent of the fermion type [TJ, [|, and this suggestion has been 
borne out by the 2-generation analysis reported above. It makes practical 
sense therefore to adopt the same position here, especially since it raises the 
stakes and makes the present analysis an even more stringent test for the 
rotation hypothesis. This means in particular that the state vector, say v T , 
of the t lepton, this being the heaviest eigenstate of the lepton mass matrix, 
could be identified again with the vector r(/z) taken at the scale /i = m T , the 
location of which can readily be determined by interpolating r(/i) between 
rrib and m c , as in Figure |3|. 

Having fixed v T , one can constrain the vector r(fi) at /i = m M by the 
condition: 



v u .r(/i)| 2 = 



m u /m t ] |v d .r(/i)| 2 = m d /m b . 



(35) 




(36) 
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with: 

v M .v T = 0, (37) 

or in other words: 

1 - |r(m At ).v r | 2 = m M /m T . (38) 

This gives 4 solutions for r/(m^), C(ra M ), corresponding to respectively the 
two signs of r(m At ).v T and the two signs of £,{m^) = ±^/l — rjim^) 2 — £(m M ) 2 . 
These 4 solution are, however, widely separated, so that only the solution 
with r(m At ).v T and £ (rn^j both positive is shown in Figure [| which, by in- 
putting the empirical values of m T = 1777 MeV and m M = 105 MeV taken 
from ||, gives as the allowed region for r(m At ) a narrow band on the \x plaque- 
tte approximately parallel to the £-axis, the width of the band representing 
the error on v T obtained from the above interpolation. For the rotation hy- 
pothesis to be valid, the trajectory for r(/x) is required to pass through this 
band at fi = m M . 

The above information on the vector r(m At ) determines also to a fair ap- 
proximation the state vector v^, the latter being constrained by the "leakage" 
mechanism to lie on the plane containing v r and r(m M ) and to be orthogonal 
to v r . That this is so can be seen as follows. As noted already, the vector v T 
being near the vector r(m c ) and therefore lying very nearly on the £ ry-plane, 
the allowed band for r(m^) defined by (|38| ) is very nearly parallel to the 
£-axis so that the second component of r(m fJi ), namely rj (m M ), is very well 
determined, as is depicted in Figure |6|. By interpolating with a curve drawn 
through the four quite accurate points for respectively t, b, c and /i, one can 
then get a fair estimate for the value of T](jn s ). Hence from Figure [5], one 
can read off the corresponding value for the third component C{m s ) which on 
insertion into Figure then allows for an extrapolation to the \i mass scale 
to give an estimate for the value of C{ m ii)i which though rough, being in any 
case small, is sufficient for our purpose. Having then obtained the vector 
r(m M ), the state vector of //, namely v M , is also determined by the conditions 
stated at the beginning of the paragraph. At the same time, of course, the 
state vector v e of e is also determined by orthogonality to both and v T . 
The actual numerical values we so obtained for the charged leptonic triad 
which we shall use later for our analysis are as follows: 

v T = (0.9975,0.0700,-0.0015), 
= (-0.0654,0.9516,0.3003), 
v e = (0.0224,-0.2995,0.9538). (39) 
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This determination of the charged lepton triad on which the analysis in the 
remainder of this section depends is about the best that one can do for the 
moment but is obviously not as accurate as one could wish. Nevertheless, 
as we shall see, it still serves its purpose in allowing us to extract some 
interesting information on r(/x) for the low /i region. 

First, according to the "leakage mechanism", the mass of the electron is 
given by: 

|r(m e ).v e | 2 = m e /m T , (40) 

which, as for (|36| ) and for the same reasons, gives 4 solutions, 2 of which 
corresponding to negative values for £ (m e ) lie outside Figure |], leaving two 
which are represented by respectively the line drawn on the e-plaquette and 
another line (not shown) nearly parallel to the first but 0.035 units lower 
(and hence almost coinciding with the one shown). Again, for the rotation 
hypothesis to be valid, the trajectory for t(ji) has to pass through one of 
these 2 line at \i — m e . In obtaining these lines, we have of course input the 
well known value of 0.51 MeV for the mass of the electron. 

Secondly, the MNS || lepton mixing matrix elements and U e s, as 



studied in oscillation experiments on respectively atmospheric neutrinos [13 



14 and reactor neutrinos such as are given by the inner products: 

= v^.va, (41) 

and: 

U e3 = v e .v 3 , (42) 

with v 3 = r(rn V3 ) being the state vector of the heaviest neutrino z/ 3 . With 
the \i state vector as determined above in (j3~9|), one obtains by inputting the 
experimental range for (f/^l 2 of about 1/3 to 2/3 fi~3| , |i~4j] again four solu- 



tions for the allowed region, three of which lie outside Figure |J leaving one 
(corresponding to and £(m U3 ) both positive) which is represented in the 
figure by the area bounded by the two near vertical lines on the z/3-plaquette. 
Similarly, with the e state vector as determined in (|39|), one obtains by in- 
putting the experimental bound \U e3 \ 2 < 0.027 ||15| , four solutions for the 
allowed region, but this time only two (corresponding to £,{rn UA ) negative) 
can be ignored being outside Figure f|, the other two being adjacent merge 
into one as represented in the figure by the area bounded by the two near 
horizontal lines on the ^3-plaquette. The consequent allowed region for the 
vector r(m V3 ) is thus represented by the roughly rectangular area shown, 
where we have put m 2 3 ~ 3 x 10~ 3 eV 2 as preferred by (13], [T3| and ]16 . 
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Finally, the mixing element U e 2 as inferred from solar neutrino experi- 
ments is given as: 

U e2 = v e .v 2 , (43) 

where v 2 is the state vector for the second heaviest neutrino z/ 2 which is by 
definition orthogonal to V3. Following thus the same procedure as in the 
preceding paragraph, one can determine the allowed region for the vector 
V2 by inputting the bounds on V3 as obtained above and the experimental 
bounds on U e2 |13], [T7|, [LB]. However, to extract r{m V2 ) from this, one would 
need m U2 which is experimentally still largely unknown so that the above 
information on v 2 cannot readily be presented in Figure £|. But, as we shall 
see, there is another way of displaying this information. 

5 Discussion 

The allowed regions for the vector r(^) for various scales fi displayed in 
Figure [| and its projections Figures |5|-|7|, represent all the information on r(/z) 
that can be extracted at present from fermion mass and mixing parameters, 
apart from the u, d masses and the solar neutrino angle ?7 e2 already noted. 
This information was extracted on the assumption that both fermion mixing 
and lower generation masses arose solely as consequences of the rotation of 
the mass matrix, under which circumstances the rotation is encapsulated 
entirely in the rotating vector r(/z), as was explained in equation (p~0|) or 
(|TT|). The definition of masses, state vectors, and mixing matrices, which 
for a rotating mass matrix is delicate, followed the prescription given in the 
introduction, which seems to us the natural one and is to our knowledge the 
only one available in the literature. Apart from these, no other theoretical 
input or assumption has been introduced. Hence if the rotation hypothesis 
set out above is correct, then the allowed regions should line up along some 
smooth 3-D curve from the heaviest t to the lowest v^. Indeed this is seen 
in the above-quoted figures to be the case. In the high energy region, say 
down to the fi lepton mass scale, where the allowed regions are mostly small, 
the alignment is seen to be quite accurate, not only in the projection of 
Figure | on to the /i?7-plane as already noted in Section 2, but also in the 
two other directions as seen in Figures |5] and [7|. Below the fi lepton mass 
scale, the allowed regions are larger and the constraints not too stringent but 
they are seen still to be thoroughly consistent with alignment on a smooth 
trajectory spanning some 13 orders of magnitude in energy. This looks to 
us nontrivial and lends direct empirical support to the rotation hypothesis 
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which is entirely mo del- independent and free from extraneous theoretical 
bias. This is the main conclusion we set out to establish. 

Next, we turn to consider the shape of the trajectory traced out by the 
data in Figure f|, to understand which more theoretical input will of course 
be needed. We shall do so with reference to our DSM scheme where the 
whole rotation idea was first proposed and in which a perturbative method 
for calculating the rotating trajectory was suggested and carried out already 
to 1-loop order [12|, |§. First, we notice from Figures ^ and [7] that r(/x) 
seems to be approaching asymptotic limits for both /i — > oo and fi — ► 0, and 
thus indicative of rotation fixed points at these scale values. Now, these fixed 
points are predicted by the DSM scheme, and they occur there by virtue only 
of the inbuilt mechanism in the model for driving the rotation and are thus 
independent of the adjustable parameters of the model. The detailed shape 
of the trajectory, however, does depend on the parameters of the model, of 
which there are 3, with 2 giving the initial direction of the trajectory and 1 
other governing the rotation speed. The calculation done already a few years 
ago [H] with the 3 parameters in the model fitted to the mass ratios m c /m t , 
m^/m T and to the Cabibbo angle is reproduced in Figures ||-|7[ Although 
this result has never been explicitly presented before, it can be inferred from 
e.g. Figure 3 of |J and transformed to the present frame ( |2"8"D from the frame 
used there through the vectors: 

v< = (0.9999,0.0117,0.0008); 

v c = (-0.0110,0.9148,0.4038); 

v u = (0.0040,-0.4038,0.9149) (44) 

obtained from the previous calculation ||. It is seen to agree very well with 
the newly extracted information down to the \x mass scale. In particular, 
the rotational fixed point predicted by the DSM at /x = oo is seen to be 
fully consistent with the data. Below the \l mass the DSM curve calculated 
to 1-loop order begins to deviate from the regions allowed by experiment. 
For example, on the e-plaquette in Figure |, the DSM curve if exact should 
hit the allowed line at \i = m e but, as indicated by the little cross, it hits 
the plaquette instead at some distance from the allowed line. This deviation 
represents the difference in the mass of the electron as predicted by the old 
calculation from its true value, i.e. 6 MeV instead of 0.51 MeV. Such a 
deviation is of course expected, since at lower scales, the vector r(^t) moves 
further and further from the high energy fixed point predicted by the scheme 
so that the 1-loop calculation for the trajectory will become less and less 
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reliable. However, the 1-loop approximate trajectory from || still hits the u 3 
plaquette inside the allowed region, in other words giving correct predictions 
for the MNS mixing elements and U e ^. This is because these elements 
depend only on the vector v 3 = r(m„ 3 ) which, as indicated in Figures [| and 
0, is already near the asymptotic value. Hence, the fact the the calculation 
agrees with data for and U e % suggests that the rotational fixed point at 
\l — is correctly predicted, although the rotational curve itself near this 
fixed point is not, by the 1-loop approximation. 

In contrast, the state vector v 2 of the second heaviest neutrino u 2 rep- 
resents the tangent vector to the trajectory near the low energy fixed point 
and cannot therefore be expected to be accurately predicted by the 1-loop 
calculation of ||. Indeed, the value predicted by || for the mixing element 
U e 2 which depends on v 2 fell outside the limits set by the solar neutrino 
experiments. In our present analysis, the information on v 2 extracted from 
the experimental limits on U e 2 can be presented as a wedge-shaped region in 
Figure |5| in which the tangent to the trajectory at the low energy fixed point 
is supposed to lie, which region is estimated with a bound \U e 2\ 2 ~ 0.33 ±0.1 
favoured by present experiments [[L3], [H]. As can be seen in the figure, the 
trajectory predicted by the DSM 1-loop calculation does not satisfy this cri- 
terion. Again, as in previous cases, there are in fact four solutions to this 
allowed region, among which we have chosen to display the one which is 
nearest to accommodating the DSM 1-loop trajectory. However, this is not 
surprising since it is already expected that the 1-loop trajectory will be un- 
reliable below the /i mass scale. In that case, it may be interesting turning 
the argument around to use the information at low scale, scanty though it is 
at present, to constrain the exact trajectory if such really exists. One sees 
then that just by deforming somewhat the 1-loop curve, one would be able 
to remove both the previously noted discrepancies in the e mass and in the 
mixing element U e2 , as indicated in Figures [|, || and [7]. 

In summary, we conclude that the existing data on fermion mass and 
mixing when appropriately interpreted do support the hypothesis of a mass 
matrix rotating with changing scales, and that the rotation trajectory indi- 
cated bears a close resemblance to that predicted earlier by the DSM scheme. 

We thank Carmen Garcia Garcia for kindly helping us with the best fit to 
the data presented in Figure and Bill Scott for advising us on the neutrino 
oscillation data. 
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